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Abstract. During post- Newtonian evolution of a compact binary, a mass ratio v 
different from 1 provides a second small parameter, which can lead to unexpected 
results. We present a statistics of supermassive black hole candidates, which enables 



us first to derive their mass distribution, then to establish a logarithmically even 
probability in v of the mass ratios at their encounter. In the mass ratio range 

£■ — \ v 6 (1/30,1/3) of supermassive black hole mergers representing 40% of all possible 

cases, the combined effect of spin-orbit precession and gravitational radiation leads to 
a spin-flip of the dominant spin during the inspiral phase of the merger. This provides 
a mechanism for explaining a large set of observations on X-shaped radio galaxies. In 
another 40% with mass ratios v 6 (1/30, 1/1000) a spin-flip never occurs, while in 

C^} \ the remaining 20% of mergers with mass ratios v £ (1/3,1) it may occur during the 

plunge. We analyze the magnitude of the spin-flip angle occurring during the inspiral 
as function of the mass ratio and original relative orientation of the spin and orbital 
angular momentum. We also derive a formula for the final spin at the end of the 
5_j ■ inspiral in this mass ratio range. 
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1. Introduction 

During galaxy mergers, following a regime of slow approach due to dynamical friction, 
eventually the central supermassive black holes (SMBHs) approach each other to a 
separation of the order of 10 3 Schwarzschild radii, when gravitational radiation takes over 
as the leading order dissipative effect. The Laser Interferometer Space Antenna (LISA, 
see [T]) is expected to detect merging binary SMBHs with masses mi + m 2 < 10 7 solar 
masses (M Q ) up to redshift z « 30. A post-Newtonian approach is well suited to describe 
their forthcoming inspiral, a regime we define in terms of the post-Newtonian (PN) 
parameter e = Gm/c 2 r pa v 2 /c 2 G (e» n = 10 _3 ,£/; ra = 10 _1 ), where r and v characterize 
the orbital separation (from the center of mass) and speed of the reduced mass particle, 
G is the gravitational constant and c the speed of light. Various corrections to the 
conservative dynamics add up to 2 PN, while the gravitational radiation results in 
dissipation of energy, angular momentum and orbital angular momentum at 2.5 PN. 

The leading order conservative correction to the Newtonian dynamics in a compact 
binary, which results in a change of the orbital plane (defined by the direction Ljv of the 
Newtonian orbital angular momentum Ltv = firxv of the reduced mass particle //) is the 
spin-orbit (SO) interaction [2], [3]. The precessional time-scale (the time during which 
the normal to the orbit Ljy undergoes a full rotation) is longer than the orbital period, 
however shorter than the characteristic time-scale of gravitational radiation (defined as 
L/L, where L is the magnitude of the total orbital angular momentum). Combined with 
the leading order gravitational radiation backreaction averaged over one quasicircular 
orbit, the SO correction provides a fair approximation to orbital dynamics, explored in 
Refs. [3], [1]. 

X-shaped radio galaxies (XRGs) exhibit two pairs of radio lobes and jets [5], [6]. 
A recent review [7] summarizes the four different models explaining XRGs: galaxy 
harbouring twin AGNs, back-flow diversion models, rapid jet reorientation models, 
finally a new jet-shell interaction model. A large subset of the observations (excepting 
cases, when the jets are aligned with the optical axes of the host ellipticals jS]) are 
well-explained by the jet reorientation model, which in turn implies a spin-flip [5], [9] of 
the dominant black hole. 

The details of how this would occur were worked out in Ref . [I] . A key element was 
the determination of the typical mass ratio at SMBH mergers by a series of estimates, 
which resulted in mass ratios v = m 2 /mi = 1/30 to v = 1/3. Because the spin scales 
with the mass squared, the second spin was neglected and only the dominant spin Si 
(with magnitude Si) kept. We summarize the consequences of this model as follows. 

a) For the typical mass ratio the dominance of L over S\ is reversed as the separation 
in the binary decreases throughout the inspiral. In the last stages of the inspiral the 
spin dominates over the orbital angular momentum S\ 3> L. 

b) The angle a between the orbital angular momentum and total angular 
momentum J (with magnitude J), also the angle (3 between the dominant spin and 
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total angular momentum evolve as: 

L 
a = — -sino;>0, (1) 

$ = - sin a < . (2) 

c) The approximate expression relating a to the post-Newtonian parameter e, mass 
ratio v and initial angle a + /3 span by the dominant spin with the orbital angular 
momentum (this angle being a constant during the inspiral) is: 

sin (a + (3) . . 

tanct^ — T-fz — . (S) 

e~ x l 2 v + cos (a + (3) v ' 

(In Eq. (41) of Ref. [I] the left hand side was given as sin 2a/ (1 + cos 2a).) 

In a criticism to the work presented in Ref. [Ij, Gopakumar recently argued that "it 
is unlikely that the spin-flip phenomenon will occur during the binary black hole inspiral 
phase" [10]. This misconception comes from mixing up the instantaneous change in the 
direction of the total angular momentum, dJ/dt = [L/Jj L — f L ■ J) J ^ with its 

averaged expression (dJ/dt) = over the precessional time-scale. The angles a and f3 
are not constants during the post-Newtonian evolution, as claimed in Ref. [10], they 
rather change as given in Eqs. (JXJ)— (EJ) J^| 

In the present paper we revisit some of the arguments of the spin-flip mechanism 
and also provide more details on it as compared to Ref [I]. In Section 2 we revisit 
the typical mass ratio argument, following a recent statistics of supermassive black hole 
candidates, resulting a newly established mass distribution. We comment on how these 
findings would affect the typical mass ratio range at SMBH encounters. In Section 3 
we analyze how the spin-flip angle depends on the mass ratio and relative orientation 
of the spin and orbital angular momentum. We also derive a formula for the final spin 
during the inspiral. Finally we present our Concluding Remarks. 

2. The sky in black holes: new statistics, consequences for the mass ratio 
at SMBH encounters and chances of the spin-flip during the inspiral 

First we summarize the arguments of Ref. [4J on the mass ratios at SMBH encounters. 
The mass distribution &bh(Mbh) of the galactic central SMBHs in the mass range 
3 x 10 6 t3x 10 9 M is well described by a power-law with an exponential cutoff, 
but for our purposes can be adequately approximated by a broken power-law [TT]-[T3] 
(confirmed by an observational survey |14|). The break is at about 10 8 M Q . In agreement 
with these arguments and observations we assume §bh{Mbh) oc Mbh-i with k £ (1,2) 
below, and ^bh(Mbh) oc Mgjj, with h > 3 above the break. Then the probability for a 
specific mass ratio arose as an integral over the black hole mass distribution, folded with 

| Only when the total and orbital angular momenta are aligned, become the angles a and /3 individually 
constant, as they identically vanish. Therefore in the aligned configuration no spin-flip could ever occur 
by the combined mechanism of SO precession and gravitational radiation. 
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the rate F to merge, and by adopting the lower values of the exponents. For the merger 
rate we assumed that it scales with the capture cross section S (the dependence on the 
relative velocity of the two galaxies was neglected, as the universe is not old enough for 
mass segregation). For the capture cross-section we assumed S oc z/ -1 / 2 , motivated by 
the following arguments: 

• for galaxies an increase with a factor of 10 in radius (10 2 in cross-section) accounts 
for an increase with a factor of 10 4 in mass (from the comparison of our Galaxy 
with dwarf spheroidals [T5]-|16]. 

• there is a well established correlation between the SMBH mass and the mass of the 
host bulge [T7] . 

• the mass of the central SMBH scales with both the spheroidal galaxy mass 
component and the total, dark matter dominated mass of a galaxy [18J. 

As a result of these considerations we have found that most likely the mass ratio is 
in the range v e (1/30 , 1/3). A typical value to consider would be v = 1/10, thus one 
of the SMBHs being 10 times as massive as the other. 

New work on the statistical analysis of 5,895 NED candidate sources [19] has been 
carried out in the mass range from 1O 5 M to above 1O 9 M . Below about 1O 6 M 
all candidates are probably compact star clusters, however the rest are likely SMBHs. 
This work shows that the SMBH mass function is a broken power law with Mg# at 
low masses, and Mg| at high masses, with a break near 1.25 x 1O 8 M ; this general 
behaviour has been long known, and has now been rederived with a very large sample. 
The key difference with respect to previous work was the careful attention paid in order 
to have equal probability for detecting a SMBH in a galaxy, regardless to the Hubble 
type. The mass distribution of the SMBHs is represented on Fig [U This particular 
distribution can be interpreted in the context of the merger model [20] with a merger 
rate scaling as (mass) +2 , very much stronger than what we favored in Ref [I]. The 
extreme mass dependence describes well a M^% black hole mass distribution consistent 
with the high end of the mass distribution; on the other hand a mass dependence of 
the merger rate close to (mass) + ^ 3 , suggested by gravitational focusing arguments [20] 



describes well the lower mass distribution nearer to Mg#. It remains to be seen, whether 
all details of the mass function can be understood using either of these mass ratio 
dependences. Of course these simple merger rate calculations assume an environment 
without cosmological expansion. However, for the densest part of the cosmos the local 
expansion is very weak [21], and that is where most of the mergers occur. 

However, for the determination of the typical mass ratio the essential result is only 
slightly changed. A merger rate running with mass" 1-2 analytically gives a M^% mass 
function (see [20]), as observed; we use this rate to estimate here the typical mass 
ratios for high BH masses. Redoing the integrals of Section 2 of Ref. [I] with k = 2, 
h = 3 (denoted there a, j3) and £ = 2 (as in [20] , so much more extreme than what was 
assumed in [1]), then all four integrals are still dominated by the lower bound; only the 
second of the integrals has q = z/ _1 in its lower bound, and so the four integrals have the 
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g-dependencies of q°, q +1 , q~ l and again q~ l . We can ignore the second integral, since 
it all refers to lower masses merging with lower masses. The most important integrals 
are those combining a SMBH above the break with a SMBH either below or above the 
break. Then the distribution in q is found as q~ l , a logarithmically even distribution in 
(dq)/q over a range of q from 1 to 1000, so a logarithmic average of 30. Weighting the 
two parts of the distribution, the larger mass ratios are favored, which would skew the 
logarithmic average of the mass ratio to q > 30, thus v < 0.03. 

The logarithmically even distribution means that the mass ratio ranges v from 1 
to 1/3, from 1/3 to 1/10, from 1/10 to 1/30, from 1/30 to 1/100, from 1/100 to 1/300 
and finally from 1/300 to 1/1000 are roughly equal likely. A glance at Table 1 of Ref. 
[I] shows, that concerning the behaviour of the ratio of the dominant spin and orbital 
angular momentum magnitudes, we have three regimes: 

(1) v G (1/3, 1) when S\ < L throughout the inspiral, 

(2) v G (1/30, 1/3) when the initial S\ < L is reversed to Si > L during the inspiral 
and 

(3) v E (1/1000, 1/30) when Si > L holds throughout the inspiral. 

For the mass ratio ranges (1) and (3) no spin-flip can occur during the inspiral, 
while for (2) it should. For (1) there is chance for a spin-flip to occur during the 
plunge, as some numerical simulations have already found this for equal masses [22J. 
For (3) by contrast there is no possibility for a spin-flip by the combined mechanism 
of SO precession and gravitational radiation. These mass ratio ranges then occur with 
(1) 20%, (2) 40% and (3) again 40% probability. This means that the spin-flip still 
typically occurs during the inspiral. 

3. Spin-flip angle distribution 

In this section we will present an analysis of the spin-flip angle occurring during the 
inspiral phase in the mass ratio range v € (1/30, 1/3) as a function of the mass ratio. 

The spin-flip model be understood as follows. Initially the galactic SMBH has 
conserved spin, along which the primary jet can form. When the two galaxies collide, the 
SO induced spin precession starts, while gravitational radiation is diminishing the orbital 
angular momentum. The direction of the total angular momentum stays unchanged. 
The constancy of J over the precessional time-scale is due to the fact, that the change 
in the total angular momentum J =LL is about the orbital angular momentum, which 
(disregarding gravitational radiation) undergoes a precessional motion about J. This 
shows that the averaged change in J is along J (simple precession, [3]). This conclusion, 
however, depends strongly on whether the precessional angular frequency Q p is larger 
than a and ft. Indeed, if these are comparable, the component perpendicular to J in the 
change J =LL will not average out during one precessional cycle, as due to the increase 
of a it can significantly differ at the beginning and at the end of the same precessional 
cycle. Such a situation would occur, when the spin and the orbital angular momentum 
are of comparable magnitude (Si ~ L, a regime through which a binary with typical 
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Figure 1. Aitoff projection in galactic coordinates of 5,895 NED SMBH candidate 
sources. The complete sample is complete in a sensitivity sense, in order to derive 
densities one needs a volume correction. In the electronic version the colour code 
is Orange, Green, Blue, Red, Black corresponding to masses above 10 5 M Q , 10 6 M©, 
10 7 M Q , 10 8 M Q , 10 9 M©, respectively. With the exception of the less numerous first 
range (Orange), representing compact star clusters, the rest are SMBHs. 



mass ratio will pass through during the inspiral) and also roughly antialigned, a low 
probability regime known as transitional precession. During simple precession Eqs. (CQ)- 
02]) governing the evolution of the angles a and /3 also hold in an average sense over the 
precessional time-scale. In what follows, we assume simple precession. 

The magnitude of the spin is unaffected by gravitational radiation, therefore by the 
simple rule of addition of vectors the spin has to align close to the J direction. The 
second jet then can start to form. In the intermediate phase when the spin precesses, 
instead of jet formation the precessing magnetic field creates a wind, sweeping away the 
base of the old jet, which in many cases can be observed. 

3.1. Spin and orbital angular momentum orientations, final spin formula 

The key equation to start with is Eq. (J3J). In order to see the validity of this equation, 
also to generalize it to the cases of non-extreme rotation, we need to evaluate 
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Here V is some characteristic rotational velocity, R the radius of the SMBH (of the order 
of its Schwarzschild radius) and Xi G (0, 1) is the dimensionless (xi = 1 for extreme 
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rotation). Therefore 

§ - Xis 1 ^' 1 • (5) 

Next we express J/ Si first from J = Lcosa + Si cos/3 by rewriting = (a + 0) — a, 
and secondly from the equality of the projections perpendicular to L of the total and 
spin angular momenta J sin a = S\sm(a + 0), so that we can equal them. By also 
employing Eq. fl5]) and basic trigonometry we obtain 

- sin (q + 0) 

tan a « — , — -jz — . (6) 

X{ e- x l 2 v + cos [a + 0) 

This is the desired generalization of Eq. (jHJ) . 

It is worth to note that when applied to the final configuration £/j n , Eq. (151) 
also stands as a formula for the final spin at the end of the inspiral, giving the polar 
angle of the final spin afi n with respect to the axis J in terms of the mass ratio, spin 
magnitude and angle span with the orbital angular momentum. Related formulae based 
on numerical runs were advanced in Refs. [23] . These results are not immediate to 
compare with ours, as Eq. (JH]) could at most be applied at the end of the inspiral; 
although in the mass ratio range where it is valid, one would intuitively expect that as 
not much orbital angular momentum is left at the end of the inspiral in comparison with 
the dominant spin, the direction of the latter will not be significantly changed during 
the plunge. 

3.1.1. Particular cases. There are three particular configurations worth to mention: 

i) The spin is aligned with the orbital angular momentum: a + = 0, thus from 
Eq. ([6]) a = and there is no room for any spin-flip. This would be the situation for 
perfectly wet mergers, which align the spin with the orbital angular momentum. 

ii) The spin is anti-aligned with the orbital angular momentum, a + = ir. 
Therefore depending on which of the Si and L are larger, the angle a is either or 

TV. 

iii) For the parameter ranges when the denominator vanishes a + = 
arccos f — Xi 1 ^ 1 ^ 2 ^), from Eq. ([6]) we obtain a = 7r/2, therefore (3 is also determined. 

3.1.2. Discussion as function of mass ratios. Keeping in mind that due to Eqs. (CO)-® 
the angle a+(3 is a constant during the inspiral (a parameter), and the dimensionless spin 
Xi behaves similarly, also regarding the mass ratio as a third parameter characterizing 
the particular merger, the angle a in general remains a function of e, thus it evolves 
together with the orbital separation r and velocity v. 

For the mass ratio v = 1/10 we have (Si/L) in w Xi e iL l>l = 0.316xi and 
(Si/L)r in ~ Xi £ fin u ~ l = 3.162%!. As tana < S\/L and tan/3 < L/S\ (the equalities 
arising when the spin and orbital angular momentum are perpendicular) we have 
tana in < 0.316xi and tan/3/ in < Q.^l^Xi 1 - F° r extreme rotation (xi — 1) we obtain 
a in , P fin < 0.316 = 18.105°. 
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For v = 1/3 we obtain tana in < (S\/L) in s=s Xi £ in V ~ X = 0.095xi and for extreme 
rotation «j n < 0.095 = 5.44°. In fact at the beginning of the inspiral this latter condition 
(meaning that the orbital angular momentum is roughly the total angular momentum) 
holds in the whole mass ratio range v G (1/3, 1). Under these conditions Eq. ((6]) can 
be approximated as 

Win ~ Xi £ ln V ~ X sin/3 in = 0.032xi^ _1 sin f3 in . (7) 

For v = 1/30 we obtain tan/3/ in < (L/S\)f in ~ X\ l£ fin v = 0.105X] -1 an d f° r 
extreme rotation (3f in < 0.105 = 6.04°. In fact at the beginning of the inspiral this latter 
condition (meaning that the dominant spin is roughly the total angular momentum) 
holds in the whole range v G (1/1000, 1/30). Under these conditions Eq. Q can be 
expanded (to first order in /?/*«,, with Xi £ fin v °f the order of p) as 

0fin ~ Xi l£ ]in 2p sin a fm = S.162xi 1 iy sma fin . (8) 

For slowly rotating SMBHs with Xi ~ 0-1 the above formula would hold only in the 
range v e (1/1000, 1/300) . 

3.2. The spin-flip angle during the inspiral 

A minimal value for the spin-flip angle a arises by forming the difference between the 
angles (3, characterizing the orientation of the spin with respect to the inertial direction 
J. Thus 

""min Pin Pfin Oifin CYin ■ \J ) 

In the second equality we have used that ai n + (3i n = afi n + /3fi n . 

However we have to take into account, that the above is only true in a 2-dimensional 
picture. In reality the 3-dimensional SO precession will complicate the situation, and 
the above angle emerges only if the number of precessions during the inspiral is an 
integer multiple of 2tt. If instead is of the type (2k + 1) 7T the spin-flip angle will be 
maximal, to be calculated as 

Omax = Pin + Pfin — llT = 2 (a in + /3 in ) - In - (a in + a fin) , (10) 

where I = if /3 in + (3f in < n and I = I if tc < /3 in + (3f in < 2ir. 

The difference between cr max and o" min is due to the fact, that the realignment of 
the spin along J is not perfect. The closer S{ m is to J, the less their difference ought 
to be due a more perfect alignement, therefore a m3iX — cr m ; n = (3fi n — (3 in should go to 
with decreasing v. 

For generic mass ratios v G (1/30,1/3), Eqs. (jSJ), (EJ) and ( fTDl give the range of 
allowed spin-flip angle for each relative orientation a + (3 of the spin with respect to the 
plane of motion and each xi- The generic numerical solution for a m i n in the case xi = 1 
is represented on Fig [2] as function of the relative orientation of the spin and orbital 
angular momentum a + (3 and mass ratio v. For a given mass ratio the spin-flip angle 
has a maximum shifted from 7r/2 towards the anti-aligned configurations. The figure 
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confirms the prediction, that significant spin-flip will occur during the inspiral in the 
mass ratio range v G (1/30, 1/3). For mass ratios smaller than 1/100 the spin does not 
flip at all, as the infalling SMBH acts as a test particle. 
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Figure 2. The spin- flip angle c m i n as function of the relative orientation of the spin 
and orbital angular momentum a + j3 (a constant during inspiral), and mass ratio v. 
For a given mass ratio the spin-flip angle has a maximum shifted from 7r/2 towards the 
anti-aligned configurations. The mass ratios v = 1; 1/3; 1/30 and 1/1000 are located 
on the log^ -1 axis at 0; 1.09; 3.40 and 6.91, respectively, confirming the prediction, 
that a significant spin-flip will occur in the mass ratio range v £ (1/30, 1/3). For mass 
ratios smaller than 1/100 the spin does not flip at all, as the infalling SMBH acts as a 
test particle. 



4. Concluding Remarks 

In light of the new data on a large sample of SMBH candidates we have established 
that the mass ratios obey an even logarithmic distribution in v. In the mass ratio 
range v G (1/30, 1/3) of SMBH mergers representing 40% of all possible cases, we have 
investigated the SO precession driven conservative and gravitational radiation driven 
dissipative contributions to the orbital evolution during the inspiral, averaged over the 
precession time-scale. In this mass range the ratio of the dominant spin magnitude and 
orbital angular momentum magnitude S/L changes from less than 1 to larger than 1 
during the inspiral. As the direction of the total angular momentum is unchanged on all 
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time-scales larger than the precession time-scale, while the magnitude of the the orbital 
angular momentum decreases due to gravitational radiation and the magnitude of the 
spin stays constant, the spin direction has to change. The spin-flip of the dominant spin 
therefore occurs during the inspiral. If jet activity is involved, X-shaped radio galaxies 
arise by this mechanism and a large set of observations on X-shaped radio galaxies could 
be explained. 

In another 40% of the mergers with mass ratios v G (1/1000, 1/30) the spin-flip 
never occurs by this mechanism, while in the remaining 20% of mergers with mass ratios 
v G (1/3, 1) it may occur during the plunge. 

SMBH mergers of equal mass to v — 1/3 are only half as likely as the mass ratios 
1/30 to 1/3, therefore the occurrence of the spin-flip can be considered typical during the 
inspiral. We analyzed the magnitude of the spin-flip angle occurring during the inspiral 
as function of the mass ratio and original relative orientation of the spin and orbital 
angular momentum and supported by numerical analysis the theoretical prediction (Fig 
E]). We also derived a formula for the final spin at the end of the inspiral in this mass 
ratio range. 

During the inspiral the following relations among the relevant time-scales hold: tilt 
/ spin-flip time-scale > inspiral time-scale 3> precession time-scale ^> orbital time-scale 
(for all mass ratios in the typical range). Interestingly enough, the spin- flip time-scale 
for a typical mass ratio of 1/10 is only about three years, while the precession time-scale 
is less then a day [I]. Thus rapidly rotating relativistic jets coming close to our line of 
sight could produce significant variability at all wavelengths years before the coalescence. 
Therefore electromagnetic counterparts / precursors to the strongest gravitational wave 
emission are also likely to occur. 
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